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Superconducting Quantum Interference Devices (SQUIDs) and other superconducting circuits are 
limited by intrinsic flux noise with spectral density l//“ with a < 1 whose origin is believed to be 
due to spin impurities. Here we present a theory of flux noise that takes into account the vectorial 
nature of the coupling of spins to superconducting wires. We present explicit numerical calculations 
of the flux noise power (spectral density integrated over all frequencies) for electron impurities and 
lattice nuclear spins under several different assumptions. The noise power is shown to be dominated 
by surface electron spins near the wire edges, with bulk lattice nuclear spins contributing ~ 5% of the 
noise power in aluminum and niobium wires. We consider the role of electron spin phase transitions, 
showing that the spin-spin correlation length (describing e.g. the average size of ferromagnetic spin 
clusters) greatly impacts the scaling of flux noise with wire geometry. Remarkably, flux noise power 
is exactly equal to zero when the spins are polarized along the flux vector direction, forming what 
we call a poloidal state. Flux noise is non-zero for other spin textures, but gets reduced in the 
presence of correlated ferromagnetic fluctuations between the top and bottom wire surfaces, where 
the flux vectors are antiparallel. This demonstrates that engineering spin textures and/or inter¬ 
surface correlation provides a method to reduce flux noise in superconducting devices. 

PACS numbers: 85.25.Dq, 05.40.-a, 85.25.Am 


I. INTRODUCTION 


Superconducting Quantum Interference Devices 
(SQUIDs) are among the most sensitive and useful 
magnetometers.^ They are able to detect magnetic helds 
as low as 10“^^ Tesla, ^ and are currently used in a 
wide variety of applications. Examples include outer 
space tests of general relativity, detection of short-circuit 
faults in microchips, as well as several applications in 
medicine, such as measuring regions of brain activity in 
magneto-encephalography. However, this high degree 
of sensitivity also causes the SQUIDs to be sensitive to 
magnetic fluctuations intrinsic to its wires and interfaces. 

The most sensitive SQUIDs show excess flux noise (in 
addition to Johnson-Nyquist white noise) of the order 
of 1 ^$o/’\/Hz at a frequency of 1 Hz [<i)o = h/{2\e\) is 
the flux quantum]. This value has not changed in or¬ 
der of magnitude since the first measurements of flux 
noise in the 1980’s;^’^ only minor improvements are ob¬ 
served in modern devices.While flux noise is suffi¬ 
ciently low for several applications, it is still considered 
a barrier for the use of the SQUID as a quantum bit 
(the flux qubit) in a superconductor-based quantum com¬ 
puter architecture.^^ Flux noise induces dephasing and 
relaxation of flux qubits, limiting their coherence times 
to less than 20 /Lts.^^ Flux noise also degrades perfor¬ 
mance of other superconducting (SC) qubits such as the 
transmon^^ and phase qubit.® The effort to reduce flux 
noise and increase qubit coherence times has been a ma¬ 
jor source of motivation for research in improving SC 
devices. 

The flux threading a SQUID at low temperatures was 
shown to follow the Curie susceptibility law (magnetic 
susceptibility cx I/T),^^ supporting the idea that the ori- 


FIG. 1: Section of a superconducting wire with a random 
distribution of spin impurities. The wire is made of super¬ 
conducting thin-film of thickness b and lateral width W, with 
b <C W. The flux produced by a spin Si located at posi¬ 
tion Ri is given by <l>i = —F{Ri) ■ Si, with flux vector F{Ri) 
parallel to the magnetic field produced by the supercurrent 
at Ri. 


gin of flux noise was due to the fluctuation of local mag¬ 
netic moments (presumably spin impurities) near the su¬ 
perconducting wire (Fig. 1).^®“^® The identity of these 
local moments remains unknown; some of the possibili¬ 
ties are illustrated in Fig. 2, which depicts the longitu¬ 
dinal cross-section of a typical Josephson junction. Can¬ 
didates for the local moments include a variety of spin 
species: Dangling-bonds,^® interface states,^® adsorbed 
molecules,^® and nuclear spins of all atoms forming the 
materials. 

Recent theory and experiment^® provided evidence 
that low temperature/low frequency SQUID flux noise 
can be explained by a spin diffusion model, with spin dif¬ 
fusion constant measured to be of the order of 10 /rm^/s. 
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Interestingly, this value is right in the range of spin dif¬ 
fusion constants measured for lattice nuclear spins due 
to their mutual dipolar interaction.^® The typical met¬ 
als that form superconducting wires (e.g. niobium, alu¬ 
minum) all have non-zero lattice nuclear spin, that are 
expected to contribute to intrinsic flux noise up to fre¬ 
quencies of the order of 10^ Hz (the value of dipolar 
coupling between nearest-neighbor nuclei in the crystal 
lattice).^' Nuclear spin noise should be present even in 
hypothetically “perfect” devices that contain no electron 
spin centers. 

Recent experiments^®’^® claimed the observation of flux 
noise at frequencies several orders of magnitude above 
the nuclear spin cut-off frequency (10"^ Hz). Therefore, 
nuclear spins alone can not explain the origin of flux noise 
in SC devices. 

Here we describe a theory for the excess flux detected 
by SC circuits in the presence of localized magnetic mo¬ 
ments, and make explicit numerical predictions for the 
flux noise power (noise spectral density integrated over all 
frequencies) originating from electron and nuclear spins 
uniformly distributed in the surface and bulk of the SC 
wires. Our main assumption is that the superconduct¬ 
ing condensate affects the value of the spin’s magnetic 
moment only through local screening described by the 
London equations. We consider the impact of the for¬ 
mation of spin clusters (spatial spin-spin correlation), 
which is typical of electron spin systems close to a phase 
transition.^®’®® This allows the assessment of the rela¬ 
tive contributions of surface impurity electrons and bulk 
lattice nuclear spins in a variety of regimes. We present 
simple analytic expressions that allow direct comparisons 
to flux noise models and experiments. 

As we shall demonstrate, flux noise power depends cru¬ 
cially on the vectorial nature of the coupling between SC 
circuit to spins. This vectorial nature can be mapped 
out by measuring flux noise as a function of the direc¬ 
tion and magnitude of an external magnetic field applied 
along the plane of the wires. 

Previous calculations of flux noise due to localized 
spin®® were based on modelling the spin as a square loop 
of side 0.1 /im, the minimum feature size allowed by 
the finite element software FastHenry.®® More recently, 
extensive numerical studies®®’®® showed that these cal¬ 
culations greatly underestimated the value of the flux 
produced by spins located at the wire surfaces, because 
they did not take account of the singular nature of the 
spin’s dipole field. Indeed, analytic expressions for the 
flux noise power that take into account the spin’s singu¬ 
lar behaviour are still absent from the literature. Below 
we obtain these expressions and argue that they can be 
applied to arbitrary circuits made of thin-film wires. We 
shall show below that our expressions lead to much larger 
values than previous ones,® but are in close agreement 
with the most recent numerical studies.®® 

The article is organized as follows: Section H intro¬ 
duces our flux-vector model, and develops a general the¬ 
ory on how flux noise depends on flux vector and the 
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FIG. 2: (color online). Longitudinal cross-section of a typi¬ 
cal Josephson junction^’®® and expected location of the spin 
species causing flux noise. The superconducting wire is made 
of niobium, the insulator of partially oxidized aluminum 
(likely amorphous), and the substrate of silicon-dioxide. Can¬ 
didates for the spin species include electron and nuclear spins. 
Amorphous interfaces are well known sources of electron spin 
centers such as dangling-bonds®® and interface states.®® We 
also depict molecules adsorbed to the surface, that were shown 
to lead to electron traps.®® All materials have non-zero lattice 
nuclear spin that contribute to flux noise even in the absence 
of defects (e.g. Nb has nuclear spin S = 9/2). 


formation of spin clusters near a phase transition. Sec¬ 
tion HI describes our numerical dipole method to cal¬ 
culate the flux vector explicitly. Section IV describes 
our numerical results and compares them to FastHenry 
calculations. Section V describes numerical results for 
the noise power due to electron spins distributed in the 
surface of the wires, with and without spin clustering. 
Section VI describes our results for spins distributed in 
the bulk of the wire, with explicit calculations for lat¬ 
tice nuclear spins in aluminum and niobium. Section VII 
discusses the implications of our results for reducing flux 
noise and presents our conclusions. 


II. FLUX-VECTOR MODEL 

A. General case 

Consider an ensemble of localized spins labelled by 
i = 1,..., N. Each spin is located at position 12.^, and is 
described by the spin-S' operator s^. This can describe 
electron or nuclear spins, e.g. single electron impurity 
centers (S = 1/2), many-electron transition metal cen¬ 
ters (S > 1/2), or nuclear spins of lattice atoms such as 
aluminum (S = 5/2) or niobium (S = 9/2). We intro¬ 
duce the notion of the flux vector Fi = F{Ri), whose 
components F/q, describe the value of the flux for a spin 
pointing along direction a = x,y,z. The total flux that 
the spin produces on the wires forming a device is written 













3 


as 

$ = ( 1 ) 

i 

Note that the Fi are real vectors with dimensions of flux 
(si is assumed dimensionless), and $ is a scalar quantum 
operator describing flux. The sign in Eq. (1) ensures 
Fi II B{Ri), the magnetic held produced by the SC’s 
current density at the spin’s location Ri (See Section III 
below). Such a flux is directly measured in a SQUID, but 
more generally will couple to any superconducting circuit 
by producing a voltage V = —d^/dt. The problem of 
flux noise in superconducting circuits is to compute the 
thermal equilibrium noise spectral density: 

/ OO 

dt e2-/‘ (SmSm ), (2) 

-OO 

where (5$(t) = $(t) — ($), with the angular brackets 
corresponding to thermal average (^) = TrlpT^}, with 
Pt = jZ the thermal equilibrium density ma¬ 

trix, R the spin Hamiltonian, and Z = Tr 
the partition function. 

We perform a spectral decomposition of Eq. (2) by for¬ 
mally diagonalizing the spin Hamiltonian, 'H\a) = Uaja): 

^«>(/) = E ■ 

a,/3 ^ 

( 3 ) 

Any mechanism that couples the spins to the lattice^® 
or to themselves^^’^®’^® leads to hnite frequency noise. 
However, evaluating Eq. (3) for a large spin system is 
a challenging task that requires a series of uncontrolled 
approximations. 

The mechanism of frequency dependence in Eq. (3) is a 
subject of current research. In Ref. 25 the frequency de¬ 
pendence of Eq. (3) was evaluated under the assumption 
that the spins are in a paramagnetic phase with spin dy¬ 
namics governed by spin diffusion, leading to l//“ noise 
with exponent 0 < a < 1.5 dependent on frequency range 
and device geometry. Reference 29 proposed a model 
based on spin-clusters and non-diffusive dynamics driven 
by spin-spin dipolar interaction. It was argued that the 
distribution of cluster sizes gives rise to a large spread of 
spin-flip times and l/jo.ss jjgise over a broad frequency 
range, independent of device geometry. 

Here we shall focus our discussion on the total flux 
noise power. 


/ OO 

SM)df 

-OO 

= E! E ({SiQ, Sj/s}) 

~{Sia){Sjl3)j , (4) 


on the particular interaction mechanism driving spin dy¬ 
namics; thus it allows a model independent comparison of 
the role of F and the contribution arising from different 
spin species. We take a continuum limit by introducing 
the magnetization or spin density 

M{r) =-'^s.,S{r - Ri), (5) 

i 

and defining the spin-spin spatial correlation function as 

^ ({M^(r),Mp(r')}} - (M„(r)) (Mp(r')} . 

( 6 ) 

In the continuum limit Eq. (4) becomes 



( 7 ) 


with dimension d = 2 (surface) and d = 3 (bulk) to be 
considered below. 

The noise power Eq. (7) depends on spin texture 
through the spatial correlation function C°‘^{r, r'). If the 
spin system is close to a phase transition, noise power will 
show strong temperature dependence because of the for¬ 
mation of spin clusters, which are described in general by 
deviations of r') from a delta function 5{r — r'). 

In addition, Eq. (7) may be temperature dependendent 
even in the absence of spin clusters, as we show below. 


B. Noise power without spatial correlation 
(independent spins) 

There are many situations where the correlation func¬ 
tion is well approximated by a delta function, and the 
state of any spin is independent of the others. These 
include high temperature, spin textures with zero cor¬ 
relation (e.g. spin glass), and when all spins are fully 
polarized along one direction. Below we consider each 
case in detail. 

In the limit of high temperature, defined by be¬ 
ing much larger than any energy scale affecting the spins, 
we may carry out a 1/T expansion of Eq. (6) and re¬ 
tain only the leading order contribution (zeroth power 
of 1/T, which is exact at T = oo). This leading or¬ 
der contribution has all different spin configurations oc¬ 
curring with equal probability, implying {sia) = 0 and 
{siaSjp) = dij 6 a/ 3 - Further assuming that the 

spins are uniformly distributed in space with density 
leads to 


^a/3 

'~'T=oo 


S{S + 1) 


ad6{r - r')6ap, 


and the high temperature noise power 


( 8 ) 


where {A, B} = AB + BA denotes the anticommutator 
of two operators. We note that Eq. (4) is independent 




S{S + l) 


T=o 


ad / d!^r\F{r)Y 


( 9 ) 
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which is formally exact at T = oo. In the absence of 
an external magnetic field, Eq. (9) is a good approxima¬ 
tion for nuclear spins down to /iK temperatures (the en¬ 
ergy scale for nuclear-nuclear dipolar interaction between 
nearest-neighbor nuclear spins). 

If the spin system is in a phase that has approximately 
zero spatial correlation {{siaSjp) « {sia){sji 3 ) for i ^ 
j), and if they are uniformly distributed in space, the 
correlation function may be approximated by 


C:Lir,r') 

r^{R,,T) 


f°‘^{r,T)ad6{r - r'), 


(lOa) 
, ( 10 b) 


leading to the uncorrelated spin noise power 
((5$)2\ f d‘'r^F„(r)r^(r,r)E^(r). (II) 

\ / unc. J ^" 


This expression is a good approximation in at least 
three cases of interest: (I) At temperatures higher than 
the spin-spin coupling J, but lower than single-spin 
anisotropy energy (which is non-zero for S > 1 / 2 ). As¬ 
suming the single-spin anisotropy is equal for all spins, 
the function f°‘^{r,T) will depend only on temperature 
and not on r, and will differ from Sap signaling the pres¬ 
ence of anisotropy. ( 2 ) At all temperatures, when the 
spin-spin coupling alternates in sign randomly, such as in 
a spin-glass. As a result, the coarse-grained spin-spin cor¬ 
relation function will average out to zero over the length 
scales of Fa{r), but the long time averages {sia) will 
remain non-zero (non-ergodicity). A space-dependent 
function may be used to model the lack of 

translation symmetry of the spin state. In the simpler 
case of translation symmetry, a spin-diffusion model with 
/«/3(r) = 5(5+I)x(r)/ (3xo) [Eq. (9) of Ref. 25] was 
proposed to describe the frequency and temperature de¬ 
pendence of the spin-glass noise. (3) When all spins are 
polarized along the same direction, which occurs at low 
T in a ferromagnetic phase, or in the presence of a large 
external magnetic field. 

To illustrate case (I), consider a model spin Hamil¬ 
tonian with easy-axis anisotropy along ey, "Hani. = 
—KJ2i s^||) with Si|| = Si ■ ey and K > 0 the anisotropy 
energy. When J ksT K the spins will be equally 
distributed in the zLS eigenstates of Siy; after a sim¬ 
ple calculation we get /H’H = (s^y) = S'^ and = 

j±2,j_2 _ = s/ 2 , with all other correlations equal 

to zero, leading to 

S^F^{r) + ^Fl{r) . 

( 12 ) 

Depending on the direction of F(r), the noise power may 
get reduced or increased in comparison to the high T 
result Eq. (9). 

More generally, the Hermitian matrix T) can be 

diagonalized to find its eigenvectors fj and eigenvalues 


(<54>) 


/ J^ksT^K 


CTd 


d'^7 


a-y = {{si ■ — {{si ■ > 0. We may establish a 

general inequality for uncorrelated noise by noting that 
the largest eigenvalue of /“^ is smaller than S{S+1) (the 
eigenvalue of sf): 



unc. 


= cTd F{r)-fy 

7 

< 3 ((<5$)')^ 


T=00 


(13) 


If 5 = 1/2 and the spin state has translation symmetry, 
a-y will be equal to 1/4 for all 7 , and the noise power is 
simply equal to the T = 00 result Eq. (9). 

Finally, we consider case (3) in detail. At low T, the 
spins may spontaneously polarize if the spin-spin inter¬ 
action is ferromagnetic (See section IIC); alternatively 
they can be polarized with a strong external magnetic 
field. Measuring flux noise in the presence of a strong 
B field is a quite challenging experiment. For supercon¬ 
ductivity to remain unnaffected, the external B field has 
to be applied along the plane of all wire segments, with 
the magnetic length remaining larger than the thin-film 
width b. Moreover, any modulation of the critical current 
due to the presence of flux perpendicular to the Joseph- 
son junctions will have to be accounted for. Despite these 
challenges, we shall show that it is worth considering this 
experiment because it would provide a mapping of F[r) 
and the measurement of spin quantum number S. 

When the Zeeman energy scale dominates, the spin 
Hamiltonian can be approximated by 'Hzeema.n = 
— Mi with /jLi = —gpigSi. Here we use /ig = gs for 
electron spins (/tb is the Bohr magneton), and fis = —MAf 
for nuclear spins is the nuclear magneton), with g the 
(/-factor {g is quite close to 2 for most electron impurities 
but can be very different than 2 for nuclei). In this regime 
correlation between different spins will be zero and the 
functions /“^(r,T) in Eq. (lOb) can be separated into 
two kinds, /y(T) and /_l(T): 

/y(B) = {{s, ■ Bf) - {{s, ■ B)f 

If I (25 + 1 )^ 

4 \ sinh^ (B) sinh^ [(25' -|- 1)B] 

MB) = M ■ BM") = M ■ BM") 

= icoth(H) |(25'-f l)coth[(2S'-f 1)5] 

— coth(5)|, (14b) 

where B is the unit vector along the magnetic field, 
B±i,B ±2 is a set of orthogonal unit vectors perpendic¬ 
ular to it, and B = gg,sB/{2kBT) measures the strength 
of the field. The equilibrium spin polarization is given by 

(s,) = -^{(25-^1) coth [(25-k 1)5]-coth(5)}.B. (15) 
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Plugging these expressions into Eq. (11) yields 

[h\Fir)f 


« ad I <fr 

' High B 


■ (16) 


Since /±(-B) > f\\{B) for all B, Eq. (16) shows that the 
noise power gets reduced upon the application of external 
B field. In the limit B ^ 1, all spins are polarized leading 
to 


(a) 

x=- 

Z‘ 

W/2 z= 0 

X= + W/2 



® 0{E 

X 


z=-b 




(5d>)") 


pol. 


s 

\\B 2 . 


r , 


^ 2 


V 

' 71 

/ d^r 

\Fir)\^- 

F{r) ■ B 

5 




(17) 

which is quite different from the high T result Eq. (9). 
Equation (17) is the exact noise power of a fully polarized 
spin system, irrespective of whether the polarization oc¬ 
curs because of application of a B field, or due to sponta¬ 
neous symmetry breaking e.g. in a ferromagnetic ground 
state. Thus, measurements of the dependence of noise as 
a function of the direction and magnitude of the in-plane 
magnetic field, or on direction of the spin polarization 
provides information on the components of the flux vec¬ 
tor F{r) and on the value of the spin quantum number 
S. 


X 


FIG. 3: (a) Transverse cross-section of the SQUID wire (See 
also Fig. 1). The current density is shown pointing along 
the j/-direction. (b) Magnitude of the superconductor current 
density as a function of x [Eq. (33)]. 


consider the Heisenberg model with easy-axis anisotropy, 
Bs-s = ~ X ^ ( J{Ri — Rj)Si ■ Sj — K ^ ) Sj|| , (19) 


As we show below, the direction of F{r) is poloidal 
along the surface of the wire (Fig. 1), implying that 
spin polarization along a single direction can never com¬ 
pletely suppress flux noise. Note how the residual noise 
arises due to the quantum nature of the spins. For exam¬ 
ple, if the spins are polarized along x, they have equal 
probability of pointing along any of the four directions 
+y, —y, +z, —z. This uncertainty gives rise to quantum 
noise power proportional to Ff^ + B? . 

Finally, consider a spin-polarized state where each spin 
located at Ri points parallel or antiparallel to Fp, we call 
this the poloidal state.The noise power produced by this 
state is obtained by evaluating the average in Eq. (4) with 
a single product state of N spins each pointing along Fp, 
the result is identical to Eq. (17) with B —> F{r), leading 
to 


1 j * 

where K > 0 and s^n = Si ■ ey, with ey the direction of 
the anisotropy axis. We take the continuum limit with 
M{r) as in Eq. (5), and expand Eq. (19) up to the second 
order in spatial derivatives: 

-\Mf+eo E 

oc—x,y,z 

--mA. ( 20 ) 

CTd " J 

This expression is valid in the long-wavelength approxi¬ 
mation, when terms of order jV^Mp are negligible. 

The exchange interaction is assumed to satisfy 




Poloidal 


= 0 . 


(18) 



rJ{r)xaXi3 = 2 ^ 


d+2 T c 
0 


CK/S 5 


( 21 ) 


Therefore, engineering the poloidal spin texture enables 
complete suppression of flux noise. We emphasize that 
this is an exact result; it happens because in the poloidal 
state every single spin quantum fluctuation has to be 
perpendicular to Fi, producing zero flux noise. 


C. Noise power with spatial correlation 
(ferromagnetic spin clusters) 

The interaction between the spins may lead to non-zero 
spatial correlation. To see the impact on flux noise power. 


with parameters and J defined by 


.2 ^ I J d‘^rJ{r)xl 
2 / d^rJ(r) ’ 

J = g / 


( 22 a) 

( 22 b) 


Parameter models the range of the exchange interac¬ 
tion, while J models its average strength. Many com¬ 
peting interactions contribute to J(r), such as direct 
exchange (always antiferromagnetic) and RKKY (ferro¬ 
magnetic when kF\Ri — Rj I ^ 17,29,30 shall focus 
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FIG. 4: (color online). Comparison of the Flux produced by each spin as a function of spin location calculated by FastHenry 
and the numerical dipole method, for a superconducting wire of penetration depth A = 0.07 ^m, thickness 6 = 0.1 fim, and 
width W = 5.2 ^m. The coordinate x runs along the lateral width of the wire, with edges at a; = ±2.6/im. 


in the case 0 < J < oo and < cx), when the sys¬ 
tem is able to transition into a ferromagnetic state with 
M II ±e||. We shall also use Eq. (20) with a space- 
dependent ey = F{r) as a toy model that we call the 
poloidal model. When K Jcrd^oi^o/by such a model 
has the poloidal state M(r) || ±F(r) as its T = 0 ground 
state. 

We use mean-field theory;^^ this is done by subtracting 
T times the single spin entropy from Eq. (20). The re¬ 
sult is that the system develops a ferromagnetic moment 
(My) 7 ^ 0 at the critical temperature 

ksTc = + 2K) , (23) 

with correlation functions given by: 


and 

= {t-t'Y 

= &\/^. (25c) 

where ksT' = [5'(5' -|- l)/3]CTd^o J mt = 1 for T > 
and ut = 2 for T < Tc- The length scales C|b'f-L describe 
the average size for spin clusters polarized along ey, and 
the direction perpendicular to it, respectively. Note Tc > 
T', so that only ^y diverges at the transition. Eor d = 2 
we get 


^wiT) 
UiT > Tc) 
UiT<Tc) 


(^11.II 


(j±a,±P(^r,r') 


keT f d<^q 
ksT r d'^q 


(24a) 




(24b) 


where T a with a = 1,2 denote the two directions per¬ 
pendicular to ey (the correlation function between || and 
T a is zero). There are two correlation length scales ^y 


c ' diU ’’-’’0 


keT 

27Tjeo 

keT 

2 ^Jd 


Kn 


Kn 


^11 / 


, (26a) 


(5a/3,(26b) 


where Ko{x) is the modified Bessel function of the 2nd 
kind. 

The noise power for surface spins {d = 2) in 

the presence of correlations is obtained by plugging 
Eqs. (26a), (26b) into Eq. (7), leading to 




T Jd^rjd^r' Fy(r)Fy(r')iiro + J^u(r') • Tj_(r')it:o 




ut\T-Tc\ 


27r^|j J d'^r\F{r)\^ 


(27) 


T=oo 
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where we presented correlated noise divided by the T = 
00 uncorrelated noise [Eq. (9)] for convenience. Notably, 
Eq. (27) depends on model parameters only through 
and Tc- Once again, the vector nature of F{r) 
plays an important role in determining noise power: Cor¬ 
related noise tends to decrease (increase) when F{r) and 
F{r') are antiparallel (parallel) for r,r' in different sur¬ 
faces of the wire (See e.g. the top and bottom wire sur¬ 
faces in Fig. 1). 

We remark that mean-field theory neglects critical fluc¬ 
tuations and is not necessarily a good approximation for 
T close to Tc- For d = 1, 2 and iF = 0 (or S' = 1/2, when 
easy-axis anisotropy is effectively zero), critical fluctua¬ 
tions reduce the actual of the Heisenberg model to 
zero.^^ For K > 0, Tc is generally non-zero at d = 2, but 
can be substantially reduced in comparison to the mean- 
field prediction Eq. (23). However, mean-field theory is 
an excellent approximation in the region T ^ Tc, when 
critical fluctuations play no role. To see this, consider 
the T ^ Tc limit of Eqs. (25a)-(26b): We get C||)?-L ~ 
f,o\/T'/T 0, and Ko(\r — r'\/^) —>• 27r^^(5(r — r'), lead¬ 

ing to (^“^(r - r') « ksTf /{Jil)5{r - r')6ai3. Plugging 
ksT' = [S(S -|- l)/3](Td^Q J we see that this is identical to 
the T = oo limit shown in Eq. (8). 

Furthermore, mean-field theory neglects quantum fluc¬ 
tuations; this is a problem for the T Tc region, where 
mean-field theory is exact only in the limit S —>■ oo. 
When J{Ri — Rj) > 0 for all i,j, the exact T = 0 state 
of the system will have all spins fully polarized along ei¬ 
ther -|-e|| or —ey (or an arbitrary direction for S = 1/2) 
leading to 


= 0, (28a) 

^^=0^^= ^CTdd(r - r')da^. (28b) 

Hence r') remains non-zero when T —>■ 0, sig¬ 

naling the presence of quantum fluctuations. Plugging 
Eqs. (28a), (28b) into Eq. (7) we get the spin polarized 
noise Eq. (17) with B = ey. While this result disagrees 
with mean-field theory, we note that the ratio of Eq. (17) 
to Eq. (9) is proportional to 1/(5' -I- 1). When S ^ oo 
this ratio goes to zero, in agreement with the T —>■ 0 limit 
ofEq. (27). 


III. EVALUATION OF THE FLUX VECTOR: 

THE NUMERICAL DIPOLE METHOD 

When an electron travels across a closed path C in¬ 
side the SC wire, it is affected by a flux equal to = 
Aflr) ■ dl, where Aflr) is the vector potential for the 
spin Si located at Ri. For an extended wire with SC 
current density Jsc(’’) and total current Jsci the corre¬ 
sponding flux can be obtained by representing the wire by 


a set of infinitesimally thin closed paths C, and summing 
over all the paths with dl = d^rJsci'f')/Isc, a weight¬ 
ing factor that represents the fraction of current flowing 
through each path: 


= 


d^rAflr) ■ Jsc{r), 


(29) 


where the integral is over the region where Jsc ^ 0; be. 
the volume of the wire. The numerical dipole method 
consists in using the spin-dipole expression for Aflr)^^ 

gp.fXQ s,x(r- Rj) 

' 471 \r-R,\^ ’ ^ ^ 

together with an analytical approximation for Jsc{i')- 
Plugging Eq. (30) into Eq. (29) we get = —F{Ri) ■ Si 
with 


F{R) 


An 

5/U//0 

47r 


d-^r 


Vr X 


,3 (r-i^) X Jscir) 


Ir-Rflsc 
Jsc{r) 


d?i 


he \r - i2| 


(31a) 
. (31b) 


Note that we choose a sign convention so that F{R) 
points along B{R), the magnetic field produced by Jsc 
at R, and that the direction of F can be found using the 
right hand rule with the thumb pointing towards Jsc 
(Fig. 1). The corresponding Zeeman energy of the spin 
system is given by "nzeeman = -he f d‘^rF(r) ■ M(r). 
We remark that this same result can be obtained from 
the flux-inductance theorem.^® 

The current density Jsc includes contributions from 
external sources plus screening currents causing Meiss¬ 
ner effect. To find Jsc one usually has to integrate the 
London equations numerically with the help of software 
packages such as FastHenry.^^ However, there exists an 
important case where Jsc is known analytically: Thin- 
film wires of width b < X, where A is the SC penetration 
depth (A = 0.05 — 0.1 yum for most superconductors), and 
wire lateral width W that is large enough to satisfy 


A 


- 

~ bW ^ ’ 


(32) 


and ^sc ^ bPj where ^sc is the superconductor coher¬ 
ence length including the electron mean free path (this 
latter condition ensures that the London equations with 
a single length scale A is a good approximation). 

Below we present analytic results in terms of powers 
of the small parameter A. Neglecting terms that are first 
order in A, the SC wire current density can be written 

as37.38 











X (jxm) 



X (fxm) 


(a)Spin pointing along x (along the wire surface) 


(b)Spin pointing along z (perpendicular to the wire surface) 


FIG. 5: (color online). Flux produced by an electron spin impurity as a function of spin location inside the SC wire. All results 
used the numerical dipole method, for a superconducting wire of penetration depth A = 0.07 /rm, thickness b = 0.1 ^m, and 
width W = 5.2 fim. The coordinate x runs along the lateral width of the wire, with edges at x = ±2.6pm. Note how the flux 
is zero in the region inside the wire {z = —b/2 and x away from the edges). 


Jsc{r) 


2Isc 


irbW 



y 



2A 


for (1 — A) < |m| < 1, 
for |m| < (1 — A), 


(33) 


where 7 = 2(2 — e^/^)/ 7 r = 0.2236 is a numerical constant, and the coordinate u = 2x(W runs along the lateral 
width of the wire, as shown in Fig. 3. The numerical dipole method consists in using Eq. (33) in Eq. (31b) to get an 
approximation for the flux vector that neglects the feedback effect of the spins on Jsc- This feedback effect shall not be 
significant when the spins are unpolarized. Later we will confirm this expectation by comparing our numerical dipole 
method to exact integration of the London equations using FastHenry, and show that non-zero spin polarization leads 
to an asymmetry on top of this solution. 

It is straightforward to integrate Eq. (31b) explicitly for the infinitely long wire with coordinate axes shown in 
Eig. 3. For a spin located at i? = {X, 0, Z), 


F,{X, Z) 
Fy{X,Z) 
Fz{X,Z) 


gfistioW 

C In 

ix-Xf + Z"^ 

Stt j 

-1 he 

[ix-X)^ + {Z + b)^\ 


0 , 


gpsMolF Jsc(m) 

47r Jsc 


arctan 


Z + b 
x-X 


— arctan 


—)1- 

x-Xj 


(34a) 

(34b) 

(34c) 


We can integrate a few particular cases analytically; neglecting terms that are first order in A: 
• Mid-surface, {X = 0,Z = 0).' 


^ _ gt^stio 

■l-7i|v^' 

^ nW 

1-7V^ . 

II 

0 



(35a) 

(35b) 


This expression is a good approximation in the top surface away from the wire edge. For electrons with g = 2, 
and assuming \/^ < 1 we get iFmid-surface| « «-4’o (»^4>o = 10“^$o). 
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• Mid-edge, (X = ^-,Z = 

Fx = 0, 
F, « - 

• Comer, (X =^-,Z = Q): 




rVw(l-7V^ 


1-7 


A/6 




W2hW ( 


1 — tanh 


F, 


1-7\/2A] 
ffMsAio 


l + 2(i)' 




rV26VF (l - 7\/^ 


1 - 


7 (A/6) 

'^■6+(Ar 


(36a) 

(36b) 

(37a) 

(37b) 


Note how Fi close to the wire edge is sensitive to the 
ratio A/6. The mid-edge case, Eq. (36b) provides the 
maximum |Fi| for a single electron spin interacting with 

a thin-film wire. When A/6 < 1, |Fedge| ~ 

However, when A/6 ^ 1 (ultra-thin wires), the edge flux 
can be reduced by as much as a factor of three. 

At the edge, the associated local field produced by 
the squid’s current on the spin is at most |Bioc| = 
Modsc/(7rV61U). For typical /sc ~ 1 mA and VbW ~ 
1 /rm we get |Bioc| ~ 4 Gauss, that is not sufficient to 
polarize electron spins even at the lowest temperatures 
achieved in the laboratory (10 mK). In SQUIDs made of 
carbon nanotubes or other nanostructures, the value of 
single-spin flux and the local field can be much larger.^® 


IV. NUMERICAL RESULTS AND 
COMPARISON TO EINITE LOOP 
APPROXIMATION/FASTHENRY 

We now present explicit numerical calculations of 
the flux vector by numerical integration of Eqs. (34a)- 
(34c). In order to validate our numerical dipole method, 
we performed comparison calculations using the finite 
loop/FastHenry method of Refs. 15, 32. We did this by 
designing a FastHenry^^ input file that included a long 
SC wire representing the SQUID, and a small square loop 
of side 0.1 /rm representing the spin. FastHenry has the 
advantage of integrating the London equations exactly. 

Figure 4(a) shows the computed <i>jf = Fx/2 (value of 
flux for spin-1/2 pointing along x) for electron spins inter¬ 
acting with a SC wire of penetration depth A = 0.07 /im, 
thickness 6 = 0.1 /rm, and lateral width W = 5.2 /rm. 
The flux is plotted as a function of spin location X (along 
the lateral width as in Fig. 3). Each curve was calcu¬ 
lated for a different spin-wire surface distance Z. Note 
how the numerical dipole calculations are right on top 
of the FastHenry results for Z < 0.25 /rm. However, 


for Z = 0 the numerical dipole results are 50% larger. 
Figure 4(b) shows the results for = Fz/2 (flux for 
spin-1/2 pointing along z, perpendicular to the wire sur¬ 
face) . Again, we see that both calculation methods agree 
for Z < 0.25 /rm. However, at Z = 0 (wire surface) we 
find that the numerical dipole method gives a flux that 
is six times larger than FastHenry. 

Figure 5 shows the numerical dipole results for electron 
spins inside the wire (negative Z). While the flux is quite 
high at the wire surface, it decreases to zero inside the 
wire (Z = —6/2 and X away from the edges). This result 
is particularly relevant for nuclear spins, as it shows that 
nuclei inside the wire give a smaller contribution to flrrx 
noise (single nuclear spin flux is ^ 10^ times smaller than 
the single electron values shown in the figure). 

It is important to note that all calculations presented 
here assumed a SC wire of infinite length. As a result, 
the flux at the wire edges (at x = ±lU/2) are of identi¬ 
cal magnitude, i.e., the flux is symmetric with respect to 
X = 0. Actual devices will show some degree of asymme¬ 
try for the current densities at the wire edges. For exam¬ 
ple, the SQUID is a closed SC wire loop, so to satisfy the 
Meissner effect it produces higher current density at the 
inside wire edge,^^ thus minimizing the value of magnetic 
field inside the wire. Modifying Eq. (33) to include this 
asymmetry would lead to the same degree of asymmetry 
in the calculated flux vectors. 

Another source of asymmetry in Jgq occurs due to 
spin polarization in applied B fields. With B point¬ 
ing along X, the superconductor will generate a current 
difference 6/sc between the top and bottom wire sur¬ 
faces, in order to screen out the magnetic moment gen¬ 
erated by polarized spins. A simple estimate is given 
by |6/sc|A6 = \fj,B\N, with N « 2 a2LW the num¬ 
ber of spins at the top-bbottom surfaces. This leads to 
|6/scl ~ 2 fj,B<^ 2 W/b ^ 10 — 1000 /xA for 6 = 0.1 /xm and 
W = I — 100 /xm. 

While these asymmetries can be significant, they do 
not seem to modify the noise power results shown in the 
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2X/W Z/b 

(a)Top surface, spin at (X, 0, 0) or at {X,0,—b) (b)Edge surface, spin at (ibiy/2,0, Z) 

FIG. 6: (color online). Flux vector squared for electron spin at the top and edge surfaces of the wire, for A = 0.07 /rm, 
b — 0.1 fim and W = 2,5,10 /rm. The coordinate X runs along the lateral width of the top surface, while Z runs along the 
side surface. Note how the flux is sharply peaked near the wire edges. 



2X/W 



Z/b 


(a)Top surface, spin at (X, 0, 0) or at (X, 0, —b) 


(b)Edge surface, spin at (±VF/2,0, Z) 


FIG. 7: (color online). Flux vector squared for electron spin at the top and edge surfaces of the wire, for A = 0.07 /rm, 
IF = 2 ^m, and different wire thickness b. Note how the edge flux is quite sensitive to b. 


next section. This occurs because ((<5$)^) is an integral 
of over all wire surfaces; since the asymmetry in¬ 
creases the current in one region and decreases it by the 
same amount in another, the asymmetry cancels out in 
computations of the noise power summed over all sur¬ 
faces. 


V. FLUX NOISE POWER EOR SURFACE SPINS 

We will now present results for flux noise power in a 
single wire with length L, width W, and thickness b. Our 
results can be applied to devices that contain more wires 
by simply adding the noise power contributed by each 


wire segment. 


A. Without spatial correlation, no spin clusters 

Here we focus our discussion in the T = oo noise power 
with S' = 1/2 [Eq. (9)]. We recall that this constitutes the 
exact noise power for S = 1/2 in the absence of spatial 
correlation and space inhomogeneity (See section II B); 
and for S > 1/2 it can be multiplied by 4S(S -I- 1) to 
yield an upper bound on the uncorrelated noise at all 
temperatures [Eq. (13)]. 

Figures 6-8 show explicit calculations of \F\'^ for elec¬ 
trons located at the top and edge wire surfaces for a 
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2X/W 

(a)Top surface, spin at {X, 0, 0) or at {X, 0, —b) 



Z/b 


(b)Edge surface, spin at {±W/2,0, Z) 


FIG. 8: (color online). Flux vector squared for electron spin at the top and edge surfaces of the wire, for IF = 10 /rm, 
b = 0.05 ^m, and different penetration depth A. Note how the edge flux is quite sensitive to A. 


range of parameters. Here we see the extent to which 
flux noise power is sharply peaked at X = ±Wl2 for all 
wire geometries. 

The noise power due to surface electrons in a single 
wire with length L is evaluated by assuming a uniform 
area density of electrons CT2 and plugging S = 1/2 in 
Eq. (9): 


(‘f’^)Top Surf. 



4.05 (5)'/* 

(72L{^m) 

w 

(n$o)^ 

(l - 7a/^) 

2 


1/4 


,(38) 


where in the second line we approximated the exact result 
by an analytic expression (good within 1% for 10“^ < 
2bjW < 1), and in the third line we assumed = ‘2g-B- 
We hnd that the top surface noise power is to a good 
approximation independent of A. The reason for this can 
be seen in Fig. 8(a), where it is shown that only 
depends on A when X is extremely close to ±IF/2 (Note 
how the peak value increases with decreasing A, but the 
peak width remains small). As a result we find that the 
A-dependent contribution to Eq. (38) is always small. 

The edge contribution is quite different. Figures 6(b)- 
8 (b) show that it is a good approximation to assume 
constant for all Z/fe, with its value given by Eq. (36b) 


squared. Hence we get: 


)Edge Surf. ~ ^ ^^edge | 




= 3.24 



x(n$o)^- 


Comparing Eq. (39) to Eq. (38) we see that in a typical 
device (satisfying W < lO^b), the edge surface contri¬ 
bution has the same order of magnitude as the top sur¬ 
face contribution. Flux noise measurements in several 
niobium SQUIDs^® were fitted to spin-diffusion theory 
to yield (72 = 5 x 10^^ cm“^, leading to noise power 
Eqs. (38) and (39) in the range of 1 — 100 see 

Fig. 12. 

Our expressions for the noise power account for the 
spin-dipole singularity fully, giving values that are qual¬ 
itatively different and numerically much larger than 
other expressions derived in the literature. For exam¬ 
ple, Eq. (6) of Ref. 6 predicted a term proportional to 
In (A) in the top surface electron noise power of a cir¬ 
cular SQUID; in contrast, our results show that the top 
noise power is roughly independent of A, with the edge 
noise strongly dependent on A/&. 

We now compare our results to the state of the art 
numerical calculations of noise power in SQUID washers 
presented in Anton et Our Eq. (9) is larger than 
Eq. (1) of Ref. 32 by a factor of S{S + l)/*^^, which 
equals 3 for S = 1/2. After multiplying the results 
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of Ref. 32 by 3, we find that our (‘^’^)-pop+bottom 
tained after adding-up 4 wire segments forming a washer) 
is 10% smaller than theirs for the case of short wires 
(L ^ 10 fim) in a wide range of parameters. For longer 
wires {L ^ 100 fJ-m), we find that our results are as much 
as 40% smaller. This quite good agreement indicates that 
corner effects are not substantial. 


B. With spatial correlation, ferromagnetic spin 
clusters 

We now present explicit numerical results of noise 
power with ferromagnetic correlations, for the case of sur¬ 
face electrons. In order to evaluate Eq. (27) we separate 
each surface integral into the 4 surfaces making the wire: 

fd^r=f d\+ [ d\+ f d\+ f d^r. 

J «/Top Bottom «/EdgeR «/EdgeL 

(40) 

The total integral J d^r J d^r' breaks up into 16 contri¬ 
butions, of which only 5 are distinct: 2 x Top-Top (be¬ 
cause Bottom-Bottom is identical to Top-Top by sym¬ 
metry), 2xEdgeR-EdgeR, 8 x Top-Edge, 2 x Top-Bottom, 
and 2xEdgeL-EdgeR. For ferromagnetic correlations the 
first 3 are always positive, while the last two are always 
negative. 

Figure 9 shows calculations of the correlated noise 
power Eq. (27) divided by the T = oo uncorrelated noise 
Eq. (9) (this ratio is denoted corr./uncorr.), as a function 
of T/Tc, for anisotropy energy K/{a 2 ^QJ) = 0.1 along 
the easy-axis ey = z (as T —)■ 0 all spins will polarize 
along ey). The correlated noise power has a sharp peak 
at r = Tj,; within mean-field theory this peak is finite, 
and has its width and height governed by the exchange 
interaction range Quite interestingly, for ^ b 
the correlated noise can be several orders of magnitude 
smaller than uncorrelated noise in all regimes, including 
the Tc <§; T <C 00 regime where mean-field theory is 
known to be accurate. 

To shed light in this feature. Fig. 10(a) shows corre¬ 
lated noise for T ^ Tc, as a function of ^y « = ^. Here 

we see that the origin of the noise reduction is the pres¬ 
ence of spin clusters (non-zero spatial correlation). Re¬ 
markably, correlated noise is always smaller than uncor¬ 
related noise for all average cluster sizes Figure 10(b) 
explains the origin of the reduction: For f > b, the neg¬ 
ative Top-Bottom contribution (with F antiparallel) is 
activated, which nearly cancels out the Top-Top conti- 
bution when ^ b. Thus, antiparallel F for surfaces 
across from each other make inter-surface ferromagnetic 
fluctuations interfere destructively with the intra-surface 
ferromagnetic fluctuations, leading to a reduction of flux 
noise. 

While uncorrelated noise scales roughly as 
correlated noise shows different behavior governed by the 
additional length scales In the “short-range order” 
regime 6 <C ^ W, correlated noise scales instead as 



FIG. 9: (color online). Effect of a ferromagnetic phase tran¬ 
sition on flux noise power, calculated using mean-field the¬ 
ory with anisotropy energy KI{a 2 fQJ) = 0.1 along easy axis 
ey = z. The ratio of correlated noise to T = oo noise is shown 
as a function of TjTc, calculated from Eq. (27). At T ~ Tc 
a sharp finite peak appears, with linewidth and peak value 
depending on ^o, the range of the exchange interaction. For 
large ^o, correlated noise can be several orders of magnitude 
smaller than the T = oo noise in the region Tc T < oo, 
suggesting a new method to reduce flux noise in SC devices. 
The noise reduction occurs due to ferromagnetic correlation 
between the top and bottom surfaces with antiparallel F. 


®ll 

corr./uncorr. at T = 0 

X 

0.44/(S'-b 1) 

y 

1.5/(S-f 1) 

Z 

1.1/(S + 1) 


TABLE I: Exact flux noise [from Eq. (17)] for the spin- 
polarized ground state for three spin polarization directions 
ey. We assumed a long wire with Jsc || y and other coor¬ 
dinates as in Fig. 3. The zero point flux fluctuation goes to 
zero only in the classical limit S —> oo. 


CT2TA/(^6); in the “long-range order” regime with ^ » IF 
it scales as cr2TAIF/(6^^). 

When T ^ Tc the spins are polarized along ey leading 
to pure quantum fluctuation noise that is larger than the 
mean-field theory prediction. At T = 0 the noise can be 
calculated exactly using Eq. (17) for B = ey leading to 
the values calculated in Table I. 

As noted in Eq. (17), flux noise is exactly zero for the 
poloidal state. We now describe the flux noise power 
as a function of T for the poloidal model, which com¬ 
prises Eq. (20) with ey(r) = F{r) [as noted above this 
model has the poloidal state as its ground state when 
K ^ JcTd^oifo/b)'^]- The resulting correlation functions 
are formally identical to Eqs. (26a) and (26b); however, 
they describe correlations between ey directions that are 
different in different surfaces (we take e_Li = y and 
e_L2 = F{r) X y). For example, in the top sur- 
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(a)Total ratio for different W (b)Individual contributions for W = 5 fim 


FIG. 10: (color online). Noise power for the ferromagnetic model with e|| = z (same as Fig. 9) in the T ^ Tc regime, (a) 
Correlated noise power divided by uncorrelated noise power (at T = oo). The results are presented as a function of average 
spin cluster size =5 divided by wire thickness b. Note how the presence of ferromagnetic spin clusters makes the 

correlated/uncorrelated ratio smaller than one in the Tc <C T < oo regime. The origin of this effect is explained in (b), 
where each contribution is labelled by a pair of surfaces; for example, —2xTop-Bottom refers to the negative of Top-Bottom 
+ Bottom-Top contribution in Eq. (27). When ^ > 6, ferromagnetic spin-spin correlation between the wire’s top and bottom 
surfaces (where F is antiparallel) is activated, reducing correlated noise power to a value below the T = oo uncorrelated noise. 
When ^ > W, ferromagnetic correlation between wire corners (where F is also antiparallel) is activated, producing further 
reduction of correlated noise. 


face is correlated to —Mz in the edge right surface, and 
to —Mr in the bottom surface. As a consequence, the 
poloidal model shows antiferromagnetic correlation be¬ 
tween spins located in surfaces across from each other, 
with ferromagnetic correlations between spins located in 
the same surface. The resulting flux noise power as a 
function of T/T^ is shown in Fig. 11. The behavior of the 
poloidal model is seen to be quite distinct from the case 
of pure ferromagnetic correlation. The noise at T « Tc is 
more than 100 times larger, and at T > Tc correlated flux 
noise power is always larger than the T = oo uncorrelated 
case. Clearly, this occurs because for ^ > 6 the antifer¬ 
romagnetic inter-surface correlated fluctuations interfere 
constructively with the intra-surface ferromagnetic fluc¬ 
tuations; a simple consequence of having antiparallel F 
between surfaces across from each other. In contrast to 
the behavior shown in Figs. 10(a,b) noise power increases 
with increasing spin correlation length When T > T,. 
and b ^ W, correlated noise scales as a 2 L\^/h^] in 
the “long-range order” regime with ^ ^ FF it saturates 
at its maximum value proportional to a 2 LX/b^. 


with the integral running over the volume of the wire, 
and tT3 the corresponding volume density for spins. We 
evaluated this expression numerically and were able to fit 
the following expression with high accuracy in the region 
10-3 < 2b/W < 1: 


($2 


'Bulk 


S{S+1). ,2/aaTWN 1^3 [l-b (^) 




(ffMsAio) 


V b 


) 


(i-7\/2A 


(42) 

Like the case of top surface noise, Eq. (42) is independent 
of A because at the edge depends on A only in a 
very small fraction of the wire volume. This dependence 
becomes negligible after volume integration. 


We present numerical results for bulk nuclear spins in 
the typical superconductors aluminum and niobium. For 
aluminum, the ^^Al isotope is 100% abundant in nature, 
has S = hl2 and g = 1.46,^° and forms a fee lattice 
with lattice parameter 4.05 A. Thus, 0-3 = 4/(4.05 A)^ = 
6.02 X 103° ^m-3, and from Eq. (42) we get 


VI. BULK NUCLEAR SPINS 


Eor bulk lattice spins inside the wire Eq. (9) becomes 
($")Buik = / d\\F{r)\^ , (41) 

^ J wire 


($2)ai = 1.8 X 10-2 


LW 

b{g,m) 


In 


i+(#r 


1-7\/2A 




(43) 

Eor niobium, the isotope is 100% abundant in 

natural samples, with S' = 9/2 and g = 1.37,'^^ and 
forms a bcc lattice with parameter 3.30 A leading to 
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T/T, 



FIG. 11: (color online). Flux noise power of the poloidal 
model [Eq. (20) with e||(r) = .F(t')] with anisotropy energy 
K/{iJ 2 ^oJ) — 100. The ratio of correlated noise to T = oo 
noise is shown as a function of T/Tc, calculated from Eq. (27). 
At T > Tc flux noise is always larger than its T = oo value, 
in contrast to the result shown for the ferromagnetic phase 
(Fig. 9). This occurs because the poloidal model has antifer¬ 
romagnetic correlation between the top and bottom surfaces. 
At T = 0 the poloidal state has zero flux noise [Eq. (18)] sim¬ 
ilar to the result obtained in mean-field theory. Engineering 
the spin system to be close to the poloidal phase provides a 
novel method to reduce flux noise in SC devices. 


0-3 = 2/(3.30 A)3 = 5.56 x lO^^ /rm-A We get 


($2)Nb = 4.2 X 10-2 


LW 

b{^m) 



(44) 


The Nb noise power is 2.3 times larger than Al. 

Figure 12 compares the T = oo surface electron and 
nuclear spin contributions for a Nb wire loop with L = 
100 ^m, 6 = 0.1 p,m, and (72 = 5 x 10^^ cm-^. It shows 
that the contribution of bulk nuclear spins to the flux 
noise power is typically 5% of the total noise, that is 
dominated by surface electrons close to the wire edges. 


VII. DISCUSSION AND CONCLUSION 

We proposed a flux-vector model of flux noise due to 
spins in superconducting devices, and performed explicit 
numerical calculations of the flux noise power produced 
by localized electron and nuclear spins. We emphasized 
the crucial difference between electron impurities and 
lattice nuclear spins. Electron impurities are typically 
concentrated in the wire surface (where flux vector F 
is maximum), and can be substantially affected by the 
formation of spin textures above or below the critical 
temperature for a spin phase transition. Nuclear spins 
are instead distributed in the bulk of the wire, and their 
noise is well described by the T = oo limit. 


FIG. 12: (color online). Gontributions to the T = oo flux 
noise power due to surface electrons and nuclear spins, for a 
niobium wire loop with length L — 100 /rm, penetration depth 
A = 0.07 pm, and thickness 6 = 0.1 pm, as a function of wire 
lateral width W. For the surface electrons, we assumed spin 
area density (T 2 = 5 x 10^^ cm-^ (as measured in Ref. 25), and 
separated the contributions into top-|-bottom surfaces and in- 
ner-|-outer side surfaces. The nuclear spin contirbution is for 
Nb; the noise for Al would be 2.3 times smaller. 


In many cases the vectorial nature of the spin-wire cou¬ 
pling plays an essential role in determining the value of 
the noise power. This includes the case of a phase tran¬ 
sition, when additional correlation length scales appear 
in the problem (describing the average size of spin clus¬ 
ters). Even in the absence of a phase transition the vec¬ 
tor nature of F plays a role, e.g. for spins with easy-axis 
anisotropy (typical of S' > 1/2 transition metal impuri¬ 
ties) . The key experiment to directly verify the vectorial 
nature of F is to measure flux noise in SQUIDs as a 
function of the magnitude and in-plane direction of an 
applied magnetic field. Fields of 0.5 Tesla at tempera¬ 
tures of 0.1 K will polarize electron spins without affect¬ 
ing nuclear spins. Table I compares values of flux noise 
power for large B applied along different directions. Fits 
using Eq. (16) plus a field independent contribution will 
allow measurement of S and disentanglement of electron 
and nuclear spin contributions. These are quite challeng¬ 
ing measurements, but we hope they will be performed 
in the future to confirm the origin of flux noise in super¬ 
conducting devices. 

Our explicit calculations of F for thin-film super¬ 
conducting wires shows that the flux at the wire 
edges is much larger than anticipated by previous 
calculations®’^^’^® because of two reasons. First, the sin¬ 
gularity associated to the spin-dipolar field enhances the 
flux at the wire edges; second, the edge surfaces of thin- 
film wires (hitherto ignored in previous calculations) con¬ 
tributes the same order of magnitude as the top -I- bot¬ 
tom surfaces. As a result, the scaling relations for flux 
as a function of wire geometry are modified. Flux noise 
can be greatly reduced by minimizing the electron spin 
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density at the wire edge region. One might be able to 
achieve this with chemical passivation of the surface^®’^® 
or by growing the wires with layer by layer deposition 
(to reduce the number of vacancies and other defects) 
instead of the usual evaporation method. 

We also presented for the first time a realistic estimate 
for the noise power contributed by nuclear spins. Nu¬ 
clear spin flux noise has been a subject of speculation for 
several years/®’^^^^^ and we can now ascertain that it ac¬ 
counts for approximately 5% of the total flux noise power 
affecting typical superconducting devices made with alu¬ 
minum or niobium. In the future, it is conceivable that 
one will be able to design SQUIDs with much lower de¬ 
fect spin density, making lattice nuclear spins the ulti¬ 
mate source of noise to be optimized. Nuclear spin noise 
spectral density ranges from 0 Hz up to ~ 10^ Hz, the 
value of nearest neighbor spin-spin dipolar interaction. 
For very thin wires (26/IF <C 1), nuclear spin noise will 
scale proportional to 6L/IF, so that further reduction 
can be achieved by reducing b,L and increasing IF. An¬ 
other alternative would be to use materials with zero lat¬ 
tice nuclear spin, such as making the superconducting 
wires with lead (Pb). Natural lead samples have 77.9% 
of zero-nuclear-spin isotopes (^®^Pb, ^®®Pb, and ^®®Pb), 
with nuclear spin 5 = 1/2 present in only 22.1% (^®^Pb). 
Hence, a dramatic reduction in nuclear spin noise is pre¬ 
dicted for natural Pb wires. Nuclear spin noise can also 
be reduced in niobium by using metastable nuclear states 
such as ^^'"Nb, the first excited state of ®®Nb (100% nat¬ 
ural abundance). ®®™Nb has 5=1/2 with a half-life 
of 16 years,providing noise reduction by a factor of 
(ix§)/(|x^) = A. 

For electron spins, a critical experimental question is 
whether and how flux noise power depends on temper¬ 
ature. The proper answer to this question should re¬ 
veal whether or not the spin system is close to a phase 
transition. So far two experiments addressed this ques¬ 
tion, reaching opposite conclusions in different tempera¬ 
ture ranges. Ref. 8 measured a change of up to 10® in 
SQUID flux noise power when the temperature changed 
from 0.1 — 4 K; in contrast. Ref. 25 made a direct mea¬ 
surement of noise power in the lower temperature range 
0.01 — 0.1 K (see their Fig. 3), and concluded that there 
was no change within experimental error bars. It should 
be noted that even in the absence of a phase transition 
the noise power may be temperature dependent. This 
can occur because of spin anisotropy [Eq. (12)] or because 
F{r) is temperature-dependent at temperatures near the 
superconducting critical temperature [when penetration 
depth A increases with temperature and Eq. (32) no 
longer holds]. The latter effect might be contributing 
to the temperature dependence observed in Ref. 8. 

Our results provide guidance on how to reduce flux 
noise in SQUIDs and superconducting qubits by chang¬ 
ing the wire geometry and/or engineering a spin phase 


transition. Our Fig. 12 confirms the original claim®d5 
that uncorrelated noise scales roughly as (72^/IF, and 
can be greatly reduced by using wider wires with larger 
IF. However, in the presence of spin-spin correlation the 
geometrical scaling of noise is quite different due to the 
introduction of the spin cluster length scales 

We make several predictions for flux noise power in the 
presence of magnetic correlation. If the spins at the top 
surface couple ferromagnetically (antiferromagnetically) 
to the spins at the bottom surface, flux noise power gets 
reduced (increased) when the correlation length scale f 
is increased. We should mention that magnetic coupling 
across thin metallic surfaces is a well established phe¬ 
nomenon in the physics of normal metals sandwiched be¬ 
tween ferromagnetic layers. For example. Ref. 43 demon¬ 
strates that the RKKY interaction can induce magnetic 
coupling that alternates back and forth between antifer¬ 
romagnetic and ferromagnetic depending on the thick¬ 
ness and electron density of the normal metal layer. Per¬ 
haps in the future we will be able to engineer SC wires 
with similar alternation between ferro and antiferro cou¬ 
plings, allowing the corresponding control over flux noise 
induced by electron spin impurities at the wire surface. 

Moreover, we predict that flux noise power is exactly 
equal to zero for the poloidal spin texture (all spins Si po¬ 
larized parallel or antiparallel to their corresponding flux 
vector Fi). This occurs because at T = 0 spin fluctuation 
is perpendicular to the spin polarization axis, which for 
the poloidal texture is perpendicular to Fi by design; as 
a result magnetic fluctuation of the poloidal state does 
not produce flux noise. How to imprint the poloidal tex¬ 
ture to a disordered spin system is a question for future 
research. Given that Fi is parallel to the magnetic field 
generated by the SC wire, it might be possible to at least 
partially polarize the spins in the poloidal direction using 
a large device current. 

In conclusion, we predict novel methods for reducing 
flux noise in superconducting devices: We demonstrated 
with explicit numerical calculations that enhancing ferro¬ 
magnetic correlation between wire surfaces, or engineer¬ 
ing a poloidal spin texture allows the reduction of the 
flux noise power due to electron spins by several orders 
of magnitude. The remaining noise power, due to lattice 
nuclear spins, can be as much as 20 times smaller than 
the noise observed in the current state of the art devices. 
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